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FINAL  REPORT  ON  THEORETICAL  STUDY 
OF  SINGLE  AND  MULTIPLE  SCATTERING  BY  CYLINDERS 


1.  INTRODUCTION  AND  SUMMARY 

In  this  report  we  present  theoretical  results  obtained  in 
the  subject  contract  in  the  areas  of  single  scattering  and 
absorption  of  e 1 e c t roaa gne t i c  radiation  by  arbitrarily  oriented 
conductive  fibers,  and  in  radiative  transfer  by  large  aggregates 
of  such  fibers. 

It  has  been  determined  that,  in  the  case  of  single  scattering 
and  absorption,  the  theory  in  its  present  fora  is  capable  of 
producing  reliable  results  over  a  very  wide  range  of  length-to- 
wavelength  ratio.  An  asymptotic  expression  is  given  for  the 
backscatter  cross  section  in  the  large  length-to-wavelength  ratio 
limit.  In  the  small  (Rayleigh)  limit,  it  is  shown  that  a  previous 
theory  by  two  of  the  authors  is  in  agreement  with  the  more 
complex  present  theory.  Two  classes  of  cases  have  been  identified 
in  which  the  present  theory  is  in  error.  In  practical  computa¬ 
tions,  these  cases  can  be  easily  avoided  without  substantial  loss 
of  information. 

The  detailed  differential  scattering  cross  sections  have 
been  appropriately  averaged  for  input  to  the  Radiative  Transfer 
computer  program.  Representative  results  are  presented  in  this 
report.  The  program  has  been  tested  in  various  limits  including 
isotropic  scattering  over  a  wide  range  of  optical  depths.  These 
results  were  found  to  be  in  excellent  agreement  with  published 
values. 


2.  ELECTROMAGNETIC  PROPERTIES  OF  CONDUCTIVE  FIBERS 


Basic  Theory. 

This  theory  is  based  upon  a  variational  procedure  first  set 
forth  by  Ta  i  for  perfectly  conducting  wires,1  and  later  extended 
to  the  case  of  finite  condnctivity  by  Cassedy  and  Fainberg,  who, 
however,  considered  only  broadside  incidence. 2  in  the  present 
theory,  scattering,  absorption,  extinction,  and  radar  cross 
sections  are  calculated  for  arbitrary  anglea  of  incidence.  The 
results  can  then  be  averaged  over  all  anglea  of  incidence  to 
obtain  results  for  a  cloud  of  randomly  oriented  particles. 


Consider  a  plane  electromagnetic  wave  incident  upon  a 
cylindrical  wire  of  finite  conductivity  at  arbitrary  angle  of 
i  nc  i  de  nc  e  ,  Oj  ,  and  arbitrary  polarisation  angle,  ,  as  shown  in 
Fig.  1.  Assuming  the  wire  is  sufficiently  thin  for  the  current 
to  be  radially  symmetric,  one  may  consider  the  current  as  a 
current  filament  I(z)  along  the  axis.  Since  the  wire  is  assumed 
to  be  thin,  only  the  component  of  the  electric  field  parallel  to 
the  axis  will  stimulate  a  response,  and  the  integral  equation  for 
the  boundary  condition  at  the  cylinder  surface  may  be  written 

„  „  jk  7.  COS  9 : 

E  sin  0:  cos\h  eJ  1  -  I(z)Z 

o  1 


Hz' ) 


d  z'  cl<£  1 


Here,  Z  is  the  skin  impedance  per  unit  length  of  the  cylinder, 
relating  Hg  to  Ez ,  and  is  given  by^ 


where 
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Figure  1 . 


Geometry  of  the  Cylindrical  fire  With  Respect  to  the 
Propagation  Vector  k,  and  the  Electric  Field  Vector, 
E.  0 1  is  the  Angle  Between  k  and  the  z  Axis;  (pit  the 
Angle  E  Makes  With  the  k-z  Plane. 


I0(g«)  and  Il(ga)  are  Modified  Beasel  functions.  It  nay  be  noted 
that  the  real  part  of  g2  nay  becone  important  in  the  visible 
region.  In  Eq.  (1) 


ri  =  ;i  /eQ  t*le  c^arac teris t ic  impedance  of  free  space,  and 


R 


-  / 


(z-z  )2  +  4a2  sin^ 
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A  stationary  homogeneons  functional  expression  can  now  be 
constructed  for  the  current  in  terns  of  the  backscattered  amplitude 
S,  giving 


For  the  current  function  we  employ  a  series  of  harmonic  terms,  i. 


I(i)  =  I, 


N  / 
E  {* 

n=l  ' 


^n^cn^ 1 ^  +  ®n^sn^*) 


n*  sn ' 


}  • 


(5) 


Here 


fCB(z)=  cos  nkz  cos  nqx  -  cos  nx  cox  nqkz 


(6) 


a  nd 


f sn^  =  *i»  nkz  sin  nqx  -  sin  nx  sin  nqkz  ,  (7) 


where  q  =  cos  and  x  -  kb,  h  being  the  cylinder  half-length. 
This  current  function  is  summed  over  the  number  of  terms  appro¬ 
priate  for  the  ratio  of  cylinder  length  to  wavelength  of  interest, 
generally  N  2  2h/X.  For  n  =  1  ,  the  current  function  becomes 
identical  to  that  used  originally  by  Tail  and  Pedersen*.  Note 
that  both  fcn  (the  even  functions)  and  f8n  (the  odd  functions) 
vanish  at  z  =  +  h,  which  is  appropriate  if  the  current  is  to  be 


zero  at  the  ends  of  the  fiber. 


When  the  current  function  ie  substituted  into  Eq.  4,  the 
following  expression  results: 
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Here  the  terainology  of  g's,  I's.  and  y '  «  follows  the  previous 
terninology  of  Tai*,  Csssedy  and  Fainberg2,  and  Pedersen^  and  is 
s  s  follows: 
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s  2ir. 
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In  the  preceding  equations,  fcn(z)  and  fsn(z)  are  given  by 
Eqs .  ( 6 )  and  (  7  )  . 

Either  by  odd  fonction  considerations  or  by  sore  elaborate 
mathematical  arguments  it  can  be  shown  that  the  integrals  of 
Eqs.  (IS)  and  (20)  are  identically  zero.  The  remaining  integrals 
have  been  evaluated  in  closed  form,  and  the  results  listed  in  an 
earlier  report.* 


After  simplification,  Eq.  (8)  becomes 

.  N  /  2  2  N 
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S  )  n  cn  cn  n  'sn  sn  .1  n  m 

n=  1  \  m=  1 


-A  ) 

cnm  cnm 


+  B  B  (7  -  X  ) 

n  m  snm  snm 


N 

E  (A  g  +  B  g  ) 
,  n  cn  n  sn 
n=  1 


By  the  variational  technique,  the  constants  An  and  Ba  may 
now  be  evaluated  by  setting 


0 


(22) 


Equation  (22)  then  yields 
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while  Equation  (23)  yields 
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Equations  (24)  and  (25)  are  non-linear  and,  therefore,  it  is 
extremely  difficult  to  solve  for  the  An's  and  Bn’s.  To  avoid  the 
non-linearity  of  Eqs.  (24)  and  (25),  the  following  procedure 
was  used.  By  multiplying  each  term  in  Eq.  (1)  by 


fsm(z)dz, 
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where  fCB  (z)  dz,  it  given  by  Eq.  (6),  we  obtain 


scra  -  W  nE  An  (V 


-  A  ) 
cnm  cnm' 


(26) 


Similarly,  by  multiplying  each  term  in  Eq.  (1)  by 


J ^  ^cm(z)^z» 

where  fsm  (z)  it  given  by  equation  (7),  we  obtain 


cos  ^ 

~*T~ 
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Since  gcm  it  real  and  gsn  it  imaginary,  while  Ycnm>  Ytnm> 
^cnm •  and  A-snm  are  complex.  Eqt.  (26)  and  (27)  can  be  rewritten 
in  the  following  set  of  four  linear  equationa. 
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If  Eq.  (5)  Accurately  represents  the  current  in  the  cylinder, 
then  for  a  given  length  cylinder,  there  should  be  a  value  of  N 
for  which  the  values  of  the  An * s  and  Bn's  converge  as  more  terms 
are  taken.  Bqs.  (28-31)  were  programmed  for  broadside  incidence, 
but  showed  no  convergence  for  even  the  smallest  values  of  cylinder 
length.  It  was,  therefore,  concluded  that  the  series  formulation 
(Eq.  5)  is  not  an  appropriate  form  for  the  current  function. 

Next,  a  single  value  of  n  was  chosen  so  that  n  =  Integer 
( kh  +  1).  Curves  for  differential  cross  section  and  back-scatter 
cross  section  were  computed  and  plotted  for  large  values  of  kh. 
These  graphs  showed  an  unrealistic  number  of  resonances  in  the 
differential  cross-section  and  the  back-scatter  graph  did  not 
agree  with  existing  data.  It  was,  therefore,  decided  to  re¬ 
examine  the  simpler  current  function  of  Tai  where  n  -  1,  using  the 
impe  dance  Eq .  ( 2 ) . 

When  n  =  1,  Eq.  (5)  becomes 


I  =  I0  [fc(z)  +  A  fs(z)J 

where 

fc(z)  =  cos  kz  cos  qx  -  cos  x  cos  qkz 


fs(z)  =  sin  kz  sin  qx  -  sin  x  sin  qkz 


By  means  of  the  variational  technique,  A  is  determined  to  be 


gs  ('c"V 
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where  gc  ,  gs,  X.c  ,  k  s  ,  y  c  and  ys  are  given  by  Eqs.  (9),  (10), 
(11),  (12),  (16)  and  (17),  respectively  with  n  =  1.  Having  deter¬ 
mined  A,  the  general  far  field  scattered  amplitude  can  be  written 
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where  p  =  cos  0,  0  being  the  scattering  angle 


By  definition,  the  differential  scattering  cross  section  is 


then  given  as 
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The  total  scattering  cross  section 
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The  extinction  cross  section  is  defined  by  the  well-known  forward 
amplitude  theorem  as 

Att  Im  [S(u  -  0.,  6.)I 


The  absorption  cross  section  is  jnst  the  ras  power  absorbed  in 
the  scatterer  divided  by  the  ras  intensity  of  the  incident  bean. 
The  ras  power  absorbed  in  the  wire  is  given  by 
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Re(Z) 


dz 


and  the  ras  intensity  of  the  0  coaponent  of  the  incident  beaa  is 
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o 

2  n 


so  that  the  absorption  cross  section  is  given  by 


a 
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2nRe(Z  ) 


dz  . 


Finally,  the  radar  cross  section  is  defined  to  be 


°RCS  =  4"  a,e‘'0i> 


k_  |S<W| 


Equations  (37-41)  are  for  polarization  in  the  k-Z  plane,  e.g. 
0  ° . 


The  results  achieved  by  using  the  above  equations  with  n  =  1 
will  be  discussed  in  the  following  section.  It  appears  at  this 
time  that  this  theory  may  be  valid  for  kh  values  much  larger  than 


RESULTS 


As  discussed  in  the  preceding  section,  the  nse  of  n  Fonrier 
series  representation  (Eqs.  5.  6,  7)  for  the  current  function 

resulted  in  the  lack  of  convergence  of  the  series.  Considerable 
effort  was  expended  in  the  analysis  and  this  result  was,  of 
course,  disappointing. 

After  having  determined  that  the  above  non- eonve r gene e 
property  exists,  it  was  decided  that  the  n  *  1  term  of  the  current 
function,  e.g.,  Eqs.  (32),  (33),  and  (34),  would  be  utilised  in 

conjunction  with  the  surface  impedance  function  of  Wait  given  by 
Eq.  (2).  These  expressions  were  substituted  in  Eq.  (4),  and  the 
variational  technique.  Eq.(22),  with  Bn  =  A  was  utilised  to 
obtain  the  differential  scattering,  orientation  averaged  extinc¬ 
tion,  scattering,  and  absorption  cross  sections,  as  well  as  the 
radar  cross  section. 

The  objective  of  this  analysis  was  to  determine  the  range  of 
parameters  over  which  the  simplified  current  function  coaid  be 
utilized.  Questions  to  be  answered  were: 

(a)  Over  what  range  of  kh  is  the  theory  useful7 

(b)  In  the  limit  of  very  large  kb,  do  the  differential 
scattering  patterns  give  reasonable  results? 

(c)  Assuming  that  the  theory  is  well  behaved  in  the 

kh  >>  1  limit,  can  one  obtain  simplified  asymptotic 
expressions  for  the  various  cross  sections? 

(d)  How  do  these  asymptotic  expressions  compare  with 
those  which  can  be  obtained  from  infinite  length 
cylinder  calculations? 


(e)  Do  the  results  agree  with  Rayleigh  theory  in  the 
appropriate  limit? 


In  the  remainder  of  this  Section,  we  present  results  which 
demonstrate  good  asymptotic  behavior  in  the  large  kh  limit,  as 
veil  as  good  quantitative  agreement  with  published  experimental 
data  and  vith  Rayleigh  theory. 


3 .1  Differential  Scattering  Cross  Section. 

In  this  sub-section  ve  present  computed  curves  of  differential 
scattering  cross  section  per  square  wavelength  (®diff/k^)  as  a 
function  of  scattering  angle  0.  It  is  important  to  note  (see 
Fig.  1)  that  the  scattering  angle  is  measured  with  respect  to 
the  cylinder  axis,  as  is  the  angle  of  incidence  0^.  In  all 
cases,  the  electric  field  is  in  the  plane  of  the  incident  k 
vector  and  the  cylinder  axis. 

In  Figs.  2  through  9,  the  incident  wavelength  is  3.14  mm 
and  the  cylinder  radius  is  3.0  micron  ( ka  =  6  x  10~3).  The 
electrical  conductivity  is  3  x  10*  mho/m.  These  parameters  are 
representative  of  graphite  fibers  illuminated  by  a  3  mm  plane  wave. 

3.1.1  Behavior  as  a  Function  of  kh. 

Figures  2,  3,  4  and  5  show  the  behavior  of  the  scattering 
pattern  for  three  values  of  kh,  where  h  =  half  length.  In 
Fig.  2,  kh  =  1.  This  particle  should,  therefore,  scatter  like  a 
simple  dipole  and  we  see  from  the  figure  that  this  is  indeed  so. 
Note  that,  even  though  =  30°,  the  scattering  displays  a  single 
broad  lobe  which  is  symmetric  and  has  a  maximum  at  90°,  which  is 
perpendicular  to  the  axis. 

In  Fig.  3,  the  length  hss  been  increased  from  1  mm  to  5 
mm,  corresponding  to  kh  =  3.  Ve  see  that,  in  this  case,  the 
scattering  pattern  is  more  complicated.  Ve  still  see  a  little 
bit  of  the  symmetric  (dipole-like)  scattering  at  90°.  This  is  in 
the  intermediate  range  of  kh.  The  pattern  is  largely  symmetric, 
with  the  major  lobe  appearing  in  the  vicinity  of  180°  -  *  150°. 
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Figure  3  . 


Note  that,  in  the  case  of  kh  >>  1,  the  scattering  amplitude  is 
constant  in  the  axisymmetric  cones  whose  included  half  angle  is 
0,  the  differential  scattering  angle  of  the  figures.  Therefore, 
the  differential  scattering  cross  section  at  0=1 80-0  j  is  the 
same  as  that  at  0  =  180°  +  0j.  But,  this  latter  scattering  angle 
is  the  direction  of  the  incident  k  vector,  and  therefore,  is  the 
forward  scattering  direction.  It  is  well  known  that,  for  large 
kh,  we  should  expect  to  find  a  major  lobe  in  the  forward  direc¬ 
tion.  Therefore,  the  lobe  at  0  s  150°  qualitatively  agrees  with 
what  we  would  intuitively  expect  for  a  moderately  large  value  of 
kh  and  smal 1  ka . 

If  we  next  increase  the  length  so  that  kh  =  10  (a  value  well 
into  the  kh  >>  1  region),  we  see  in  Fig.  4  that  the  major  lobe 
of  Fig.  3  develops  into  a  much  more  pronounced  peak  in  the 
vicinity  of  0  ~  150°.  This  shows  that  the  forward  (0=  180°  +  0|) 
and  specular  (0=  180°  -  0^)  scattering  patterns  are  fairly  well 
developed. 

Figure  5  is  a  composite  of  Figs.  2,  3,  and  4,  with  all 
data  shown  in  correct  numerical  scale.  It  is  interesting  to  note 
that  the  peak  for  kh  =  10  is  closer  to  150®  than  that  for  kh  =  5, 
which  is  in  agreement  with  our  expectations. 

3.1.2  Behavior  as  a  Function  of  0  $ . 

In  the  next  set  of  Figs.  (6  through  9),  we  fix  kh  at  a 
value  of  kh  =  5  and  choose  three  angles  of  incidence:  0j  =  30°, 
60®.  and  90°. 

The  scattering  patten  cf  Fig.  6,  for  which  0j  =  30°,  is 
that  of  Fig.  3  and  is  repeated  for  continuity.  We  see  in  Fig.  7 
that,  when  0j  is  increased  to  60°,  the  major  lobe  shifts  to  0  ~ 
180c  -  60°  =  120°,  which  is  just  what  we  should  expect. 
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Differential  Cross  Sec t i on/Lambda^  for  kh 


At  =  90°  (Fig.  8),  ve  have  broadside  incidence  and  the 
major  lobe  has  shifted  to  90°  as  expected.  The  pattern  is 
symmetrical  about  90°  as  it  should  be. 

A  composite  of  the  preceding  three  curves  is  shown  in 
Fig.  9,  which  demonstrates  how  significantly  the  magnitudes  of 
the  scattering  cross  sections  differ  as  a  function  of  angle  of 
incidence. 

3 . 2  Backscattcr  Cross  Sections . 


We  have  made  a  number  of  comparisons  between  the  results  of 
the  present  theory  and  available  measurements.  The  data  presented 
in  this  section  are  limited  to  moderate  values  of  kh  <.  10.  It  is 
unfortunate  that  more  experimental  data  are  not  available.  In 
particular,  we  have  found  no  expet isental  data  for  lossy  cylin¬ 
ders.  Available  experimental  data  with  which  we  can  compare 
appears  to  be  limited  to  the  backscatter  cross  sections  of 
perfectly  conducting,  but  thin  (ka  <  1)  wires.  Notwithstanding, 
considerable  insight  can  be  gained  in  the  foregoing  comparisons. 


3.2.1 


ion  vs 


In  this  sub-section  we  compare  the  results  of  the  present 
theory  with  experimental  data  taken  at  Lockheed  Georgia^  for 
tungsten  wires  in  the  resonance  region.  The  measurement  frequency 
was  9.375  GHz.  Wire  diameters  ranged  from  1  to  3  mils.  The  skin 
depth  of  tungsten  is  0.04  mils  at  this  frequency  and  we  can, 
therefore,  consider  the  conductivity  to  be  infinite  for  Figs.  10 
through  16.  In  these  Figures,  the  solid  curves  are  our  theoretical 
results  and  the  dotted  curves  are  experimental  data.  These 
Figures  are  plots  of  backscatter  cross  section  per  square  wave¬ 
length  vs.  aspect  angle  with  the  E  vector  and  the  cylinder 
axis  in  the  rotational  plane.  Note  that  90°  and  270°  represent 
broadside  incidence,  while  0°  and  180°  represent  end-on  incidence. 
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Figure  10  corresponds  to  the  onset  of  the  first  ( ~  half 


wavelength)  resonance.  Note  that  the  angular  symmetry  of  the 
data  is  off  by  10  to  15  degrees.  If  this  were  corrected,  reason¬ 
ably  good  quantitative  agreement  would  be  observed. 

Figure  11  corresponds  to  the  first  resonance  (Z/X  =  0.480) 
and  we  find  excellent  agreement  between  theory  and  experiment. 
Figure  12  (  Z/X  -  0.496)  shows  good  experimental  symmetry  and 
agreement  within  10%.  In.  Fig.  13.  in  which  Z/X  =  0.525.  the 
agreement  is  better  than  10%  when  corrections  are  made  for 
experimental  assymetry. 

Note  in  Fig.  14,  for  which  Z/X  «=  0.854  the  broadside  (90°) 
peak  observed  in  Figs.  11  through  13,  has  split  and  the  RCS  is 
reduced.  Although  the  features  of  the  theory  vs.  experiaient 
agree,  the  quantitative  agreement  is  only  within  ~25%  to  30%  in 
Fig.  14  . 

Figures  15  and  16  correspond  to  Z/X  =  0.929  and  Z/X  =  1.051, 
respectively.  In  Fig.  15,  the  experimental  amplitude  assymetry 
is  about  12%,  while  the  agreement  between  theory  and  experiment 
are  within  about  18%.  In  Fig.  16,  a  new  maximum  has  formed  at 
90°  and  270°.  The  entire  pattern  is  reduced  and  only  moderate 
(20  to  25%)  theoretical  experimental  agreement  is  observed. 

3.2.2  Backscattcr  Cross  Section  vs.  kh . 

In  this  sub-section,  we  first  compare  the  results  of  the 
present  theory  with  some  e  a  x 1 y  experimental  measurements  on 
highly  conducting  wires  by  As  and  Scl.ritt®*7  with  the  results  of 
the  present  theory.  Figure  17,  the  experimentally  observed 
backscatter  cross  section  per  square  wavelength  is  plotted  as  a 
function  of  kh  for  three  values  of  ka.  Note  that  an  error  exists 
in  this  Figure.  The  solid  line  (ka  =  0.132)  and  the  dotted  line 
(ka  =  0.026)  were  obviously  inadvertently  interchanged  in  the 
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Figure  16.  Ruder  Cross  Sec t i on/Lambde^  of  u  Tungsten  Wire  us  s  Function  of 
Angle  of  Incidence,  ©i.  for  i/X.  =  1.051  end  i/d  =  884 


original  paper  by  Aa  and  Schmitt  and  the  error  was  reprodnced  in 
King  and  Vn. 

Fignre  18  ehows  the  reanltt  of  the  present  theory  for  the 
same  parameters  as  those  of  Fig.  17.  Note  that,  except  for  the 
deep  minima  of  Fig.  18,  excellent  quantitative  agreement  is 
demonstrated.  The  existence  of  these  minima  will  be  discussed  in 
Section  3.5. 

Figure  19  shows  the  behavior  of  the  backscatter  cross 
section  vs.  kh  for  four  electrical  conductivity  values  ranging 
from  100  mho/m  to  infinity.  Note  that  (a)  as  we  would  expect, 
the  backscatter  cross  sections  decrease  with  decreasing  conduc¬ 
tivity,  (b)  the  plateaus  degenerate  to  simple  maxima,  and  (c)  the 
three  very  deep  minima  persist  even  for  small  values  of  electrical 
conductivity. 

3 .3  Very  Larae  Values  of  kh . 

We  next  investigate  the  behavior  of  the  theory  when  kh 
becomes  very  large.  The  purposes  of  this  are  (a)  to  test  the 
stability  of  our  solutions  in  the  very  large  kh  range,  (b)  to 
determine,  if  possible,  asymptotic  expressions  for  the  cross 
sections,  and  (c)  to  see  if  the  solutions  appear  to  be  reasonable 
on  physical  grounds. 

In  Figures  20  and  21,  which  are  analogous  to  Fig.  17,  we 
have  plotted  (linear  plot)  backscatter  cross  section  per  square 
wavelength  vs.  kh  over  the  range  0.25  £  kh  1  100.  In  Fig.  20,  ka 
=  10-4  and  in  Fig>  21,  ka  =  10"3.  we  find  that 

(a)  the  backscatter  cross  section  appears  to  be  well 
behaved  in  the  very  large  kh  regime 


(b)  The  minima  of  Fig.  17  persist  to  at  least  kh  =  100. 


An  analysis  of  Figs.  20  and  21  gives  the  following  asymptotic 
relationship: 


Another  test  of  the  theory  is  to  observe  the  behavior  of  the 
differential  scattering  cross  section  in  the  very  large  kh 
limit.  The  results  of  these  computations  are  given  below. 

Figures  22.  23,  24  and  25  show  the  behavior  of  the  differen¬ 
tial  scattering  cross  section  per  square  wavelength  for  various 
values  of  kh.  The  cylinder  orientation  is  broadside,  with  E 
parallel  to  the  cylinder  axis.  These  values  are  kh  **  10,  15, 
25  and  100.  Note  that,  as  kh  is  increased,  (a)  the  width  of  the 
central  maximum  decreases,  (b)  the  lobe  structure  becomes  compres¬ 
sed,  and  (c)  unexpected  minor  lobes  appear  near  0°  and  180°. 

Except  for  (c)  above,  the  curves  of  Figs.  22  through  25 
appear  very  normal  and  well  behaved.  Also,  we  have  shown  that 
the  main  lobe  structure  agrees  exactly  with  the  usual  (  s  i  n2kh/ 2kh )  * 
representation. 

Figures  26  and  27  represent  an  interesting  and  unantici¬ 
pated  result.  First,  note  that  the  angle  of  incidence  is  40<>  (or 
50°  off  broadside).  In  Fig.  26  the  total  length  (2h)  of  the 
cylinder  is  an  odd  multiple  of  a  quarter  wavelength.  Specifically, 
the  total  length,  2h,  is  given  by  103.5  (1/2).  Although 
points  are  computed  for  only  every  degree  of  scattering  angle,  it 
can  be  seen  that  the  scattering  pattern  is  well  behaved,  with  a 
large  peak  in  the  forward  direction*. 

*As  was  pointed  out  earlier  in  this  report,  the  scattering  is  the 
same  at  the  angle  180°  +  0*  as  at  180°  -  0j.  Therefore,  the 
large  peak  at  180°  -  0*  can  be  interpreted  as  forward  scattering 
as  well  as  scattering  at  180°  -  0i . 
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Figure  22.  Dif ferentiel  Croce  Sec t i on/Laabda2  for  an  Infinitely 
Conducting  Cylinder,  Broadeide  Incidence  and  kh  =  10 
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Figure  23.  Differential  Cro*»  Sec t lon/Lanbda2  for  an  Infinitely 
Conducting  Cylinder,  Broadside  Incidence  and  kh  “  15 
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Figure  24.  Differential  Cross  Sec t i on/Laabda2  for  in  Infinitely 
Conducting  Cylinder,  Broadside  Incidence  and  kb  ■  25 
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Figure  26.  Differential  Croee  Section/Lambda2  for  an  Infinitely 
Conducting  Cylinder  with  =  40°  and  kh  *  163 
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Figure  27.  Differential  Croce  See t ion/Lanbda*  f0r  an  Infinitely 
Conducting  Cylinder  with  9*  *  40°  and  kh  *  162 


This  behavior  of  the  scattering  pattern  it  typical  of  nearly  all 
of  the  kh  >>  1  scattering  patterns  which  we  have  plotted.  Ve 
have,  however,  found  an  exception  to  this  behavior.  It  occurs 
within  a  very  narrow  range  of  kh  when  the  total  length  is 
approximately  equal  to  any  (even  or  odd)  multiple  of  a  half 
wavelength.  When  such  multiples  occur,  i.e.,  &  =  2h  =  NX/2  when  N 
>>  1,  we  find  a  large  peak  (equal  to  the  forward  scattering  peak) 
in  precisely  the  backsca t te r  ing  direction  0  =  0  ^ . 

This  behavior  occurs  for  all  angles  of  incidence  for  which  we 
have  made  comput a t i on s .  An  example  of  this  is  shown  in  Fig.  28. 
This  anomalous  behavior  occurs  for  less  than  5%  of  the  cylinder 
lengths  if  we  continuously  vary  cylinder  length  (0j  90°  and  kh 

>>>  1)  and  observe  the  scattering  in  the  backs c a 1 1 e r ing  (0  =  0 ^ ) 

direction.  A  further  discussion  of  this  anomaly  will  be  given  in 
Section  4. 


SecUPJf  yh 

In  the  case  of  transmission  of  electromagnetic  radiation 
into  clouds  of  randomly  oriented  conductive  fibers,  one  is 
interested  in  the  orientation  averaged  values  of  the  extinction, 
absorption,  and  scattering  cross  sections.  In  particular,  the 
cross  section  per  unit  volume  of  the  particle  is  of  importance. 
We  have  developed  a  computer  program  for  the  calculation  of  the 
above  cross  sections  as  a  function  of  the  wavelength  of  the 
incident  radiation.  The  fixed  parameters  in  the  calculation  are 
fiber  radius,  length,  and  conductivity.  One  hundred  and  twelve 
such  sets  of  plots  have  been  produced  using  the  impedance  of 
Eq.  (2)  and  the  current  function  given  in  Eqs.  (32),  (33)  and 
(34).  These  plots  are  being  forwarded  to  Dr.  Jay  Embury  at  C.S.L. 

As  examples  of  the  above  mentioned  plots.  Figs.  28,  29,  and 
30  predict  the  spectral  characteristics  of  thin  graphite  fibers 
in  the  spectral  range  0.1  micron  <  X  <  100  microns.  The  fixed 
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Figure  28.  Orientation  Averaged  Electromagnetic  Croaa  Section/ 
Lambda^  Veraua  Wavelength  for  a  Cylinder  1  x  10”^ 
Meters  in  Length 
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Lambda^  Versus  Vavelength  for  a  Cylinder  3.33  z  10  ® 
Meters  in  Length 
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parameters  are  radius  =  0.1  p ,  conductivity  =  3.3  z  1 0 l'  nho/n, 
and  total  length  =  1  p,  3.3  p,  and  10  p  for  Figs.  28,  29  ,  and  30, 
respectively.  The  dashed  lines  correspond  to  scattering  and 
absorption  in  the  Rayleigh  limit.  Note  that  the  various  electro¬ 
magnetic  cross  sections  are  averaged  over  particle  orientation. 
It  should  also  be  noted  that  the  orientation  averaged  scattering 
and  absorption  cross  sections  are  computed  by  averaging  these 
over  all  aspect  angles  and  (in  the  case  of  scattering)  all 
scattering  angles,  as  shown  by  Eqs.  (38)  and  (40).  However,  the 
extinction  cross  sections  are  computed  using  the  Forward  Amplitude 
theorem  Eq.  (39).  At  any  wavelength,  the  sum  of  the  scattering  and 
absorption  cross  sections  should,  of  course,  be  equal  to  the 
extinction  cross  section.  Since  the  extinction  cross  section  is 
computed  by  independent  means,  a  good  check  on  the  theory  is  to 
determine  whether  or  not  the  equality  5abs  +  ®sca  =  ®ext  prevails. 
Indeed,  it  can  be  seen  by  inspection  of  the  three  figures  that  the 
equality  does  exist  throughout  all  but  the  shortest  wavelength 
range  of  these. 

The  purpose  of  including  three  figures  is  to  demonstrate  the 
behavior  of  the  cross  sections  as  the  fiber  length  is  increased 
fro.ni  one  micron  to  ten  microns.  In  Fig.  28,  we  see  that  the 
absorption  cross  section  becomes  maximum  at  about  1  =  10  |i  and  is 
more  than  two  orders  of  magnitude  higher  than  the  scattering 
cross  section  for  wavelengths  longer  than  10  p. 

In  Fig.  29,  the  length  has  been  increased  from  1  p  to  3.3  p 
and  we  see  that  the  absorption  reaches  its  maximum  at  A  =  20  to 
30  microns  and  the  scattering  cross  section  is  greater  than  that 
of  Fig.  28. 

Figure  30  shows  the  behavior  when  the  total  length  is 
increased  to  10  microns.  Here,  we  see  that  the  absorption  (and 
extinction)  are  fairly  constant  over  the  range  1  p  £  A  £  100  p. 


The  above  Figure*  show  that,  increasing  the  length  of  the 
fiber  from  1  p  to  10  p  has  improved  the  absorptive  (and  extinction) 
properties  of  the  fiber  at  the  longer  infrared  wavelengths.  It 
is  also  seen  that,  as  fiber  length  is  increased,  the  ratio  of 
scattering  to  absorption  cross  sections  increases  (as  we  should 
expect ) . 

It  is  important  to  note  that  the  numerical  integrations  for 
0 s c a  and  oabs  were  carried  out  in  5°  increments  in  both  G  and 
G £ .  This  rather  coarse  increment  was  necessary  because  of 
computer  memory  limitations.  As  kh  becomes  larger,  the  scattering 
patterns  become  more  peaked,  and  the  finer  the  angular  integration 
increment  should  be.  The  use  of  a  5°  increment  is  probably  the 
reason  for  the  inequality  of  5SC#  +  ?abs  =  aCxt*  This  tentative 
conclusion  is  reinforced  by  noting  that  the  exact  specular  peak 
at  G  =  n  -  Gj  is  always  included  in  the  angular  integration  of 
the  differential  scattering  cross  sections,  and  it  is  therefore, 
weighted  more  strongly  as  the  width  of  the  peak  approaches  the 
angular  integration  increment.  We  plan  to  compute  a  limited 
number  of  plots  similar  to  those  of  Figs.  (28),  (29),  and  (30) 
using  1°  increments,  and  to  compare  the  short  wavelength  behavior 
with  that  of  Figures  28,  29,  and  30. 

Finally,  it  is  interesting  to  note  that,  in  the  long  wave¬ 
length  limit,  previous  results&,9  based  on  Rayleigh  scattering  and 
absorption  (dashed  line)  are  in  close  agreement  with  those  of  the 
present  theory . 

3 . 4  Discussion . 

The  primary  purpose  of  our  theoretical  work  on  scattering 
and  absorption  by  cylinders  has  been  to  include  ohmic  losses  in 
the  historic  formalism  developed  by  Tai*,  Cassedy  and  Fainberg^, 
and  others.  The  reason  for  this  has  been  to  provide  a  useful 
procedure  for  the  computation  of  scattering,  absorption,  and 
extinction  cross  sections  for  lossy  fibers.  It  turns  out  that 


essentially  all  of  the  available  experimental  data  with  which  we 
can  compare  the  theoretical  results  involves  wires  having  a 
conductivity  so  high  that,  for  all  practical  purposes,  it  is 
infinite.  This  is  why  a  large  portion  of  the  data  presented  in 
this  report  involve  infinitely  conducting  fibers. 

The  details  of  the  various  scattering  patterns  which  we  have 
investigated  obviously  depend  to  a  very  high  degree  upon  the 
current  function.  The  absorptive  properties,  of  course,  also 
depend  upon  the  current  function.  Ve  believe  that  good  behavior 
of  the  scattering  patterns  in  a  given  regime  justifies  the  use  of 
the  present  current  function  in  that  regime  for  the  calculation  of 
absorption  as  well  as  scattering.  To  say  this  in  another  way, 
a  careful  analysis  of  the  scattering  predictions  of  the  theory  is 
a  powerful  tool  in  determining  where  the  use  of  the  present 
theory  is  acceptable  and  where  it  is  not.  In  addition,  as 
mentioned  in  Section  3.2.3,  the  conservation  of  energy  criterion 
(oSca  +  ®abs  =  ®ext)  »l«o  serves  as  an  important  validation  tool. 

With  regard  to  the  data  which  we  have  included  in  the 
preceding  pages,  we  have  the  following  specific  comments: 

(1)  The  features  of  the  differential  scattering  data  of 
Sections  3.1.1  and  3.1.2  (Figs.  2  through  9)  appear  to 
be  as  one  would  expect,  with  the  forward  scattering 
lobe  becoming  more  pronounced  as  kh — “10.  A  somewhat 
disturbing  feature  of  these  data  is  the  presence  of 
rather  larger  lobes  near  0°  and  180°  (end-on)  than  we 
expected  to  find.  We  have  found,  however,  that,  as  ka 
is  made  extremely  small,  these  lobes  become  smaller. 

(2)  The  data  on  backscatter  cross  sections  is  very  instruc¬ 
tive.  Our  predictions  agree  very  well  with  the  Lockheed 
measurements  over  the  (rather  small)  range  of  0.480  £ 
£/X  1  1.051  in  Section  3.2.1  (Figs.  10  through  16). 


The  quantitative  and  qualitative  details  of  the  theory 
vs.  aeasnreaeit s  are  in  good  agreement. 


In  Section  3.2.1  very  good  agreement  is  demonstrated 
between  theory  and  experiment  for  the  three  ka  values 
of  Figs.  17  and  18.  However,  note  that  the  deep  minima 
appearing  in  Fig.  18  are  definitely  incorrect.  These 
zeros  in  the  backscatter  cross  section  occur  for 
arbitrarily  small  ka  values,  and  therefore,  signify 
zero  values  of  S(0),  the  forward  scattering  amplitude. 
This,  in  turn,  forces  the  conclusion  that  the  total 
cross  section  is  zero  at  these  values  of  kh,  which 
cannot  be. 

The  above  problem  can  be  directly  traced  to  the  current 
function  of  Eqs .  (32),  (33),  and  (34).  If  we  integrate 
the  current  function  over  the  length  of  the  wire  and 
set  the  integral  equal  to  zero  (e.g.  net  current  «*  0), 
we  find  that  the  relation  tan(kh)  “  kh  obtains.  Ve 
have  done  parametric  plots  of  this  equation  and  find 
that  the  zeros  of  Fig.  18  occur  precisely  at  the  values 
of  kh  which  satisfy  the  above  transcendental  equation. 
We  should  note  that,  in  re-reading  Tai's  paper,  1  we 
found  that  he  had  discovered  the  same  problem  and 
produced  the  same  transcendental  equation.  This 
is  a  definite  limitation  of  the  theory.  Aside  from 
this,  we  seem  to  be  in  good  agreement  with  experiment. 
Therefore,  taking  this  reservation  into  account,  the 
theoretical  backscatter  data  for  various  finite  values 
of  conductivity  (Fig.  19)  have  high  credibility. 

In  the  esse  of  very  large  values  of  kh,  we  see  that  the 
minima  mentioned  above  persist  out  to  at  least  kh  * 
100  (Section  3.3,  Figs.  20  and  21),  and  predictions  in 
the  near  vicinity  of  tan  (kh)  =  kh  should  be  disre¬ 
garded.  Aside  from  this,  we  find  that  the  theory  is 


well  behaved  for  extreaely  large  values  of  kb  for 
broadside  scattering.  An  asyaptotic  expression  for 
the  broadside  baekscatter  cross  section  vs  kh  was 
derived  and  is  given  by  Eq.  (35). 

The  calculations  of  differential  scattering  cross 
section  vs  scattering  angle  plotted  in  Pigs.  26  and  27 
are  very  well  behaved  and  produce  exactly  the  narrow 
forward  scattering  lobe  which  we  expected.  However,  we 
have  determined  (as  discussed  in  Section  3.3)  that, 
within  a  very  narrow  range  of  2,/A,  a  large  baekscatter 
peak  occurs  at  all  values  of  2/A  =  N,  where  N  is  any 
large  integer.  Since  this  behavior  is  essentially 
independent  of  the  angle  of  incidence,  we  conclude  that 
it  is  incorrect  and  reflects  a  limitation  of  the  simple 
current  function  which  we  are  using. 

(4)  The  computations  plotted  in  Figs.  32,  33,  and  34 

(Section  3.4)  provide  an  example  of  the  use  of  the 
theory  which  we  have  developed.  In  order  to  obtain 
accurate  results  in  the  larger  kh  range,  we  must 

obviously  use  finer  angular  increments  in  the  integra¬ 
tions  to  obtain  the  total  scattering  cross  sections. 
We  have  realized  this  for  some  time  and  are  taking 

steps  to  obtain  higher  computational  speed  and  more 
needed  computer  memory.  It  is  particularly  interesting 
to  note  the  good  agreement  in  the  Rayleigh  region 

between  the  early  calculations  of  Pedersen  and  the 
present  theory. 

It  is  recognized  that  realistic  calculations  in  the 

infrared  and  visible  portions  of  the  spectrum  require 
the  use  of  established  optical  constants  (complex 
refractive  index)  and  appropriate  modification  thereof 
in  the  very  thin  limit  of  particle  radius. 
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In  conclusion,  ve  have  found  surprisingly  good  agreement  between 
the  theoretical  predictions  and  available  experimental  data.  The 
calculations  are  well  behaved  out  to  very  large  values  of  kh 
where  we  find  no  onset  of  instability.  Two  limitations  of  the 
theory  were  identified.  These  occupy  only  a  very  small  fraction 
of  the  kh  values  in  any  given  kh  range  and  such  kh  values  can  be 
avoided  in  future  computations  without  serious  limitations  of  the 
overall  information  which  one  can  obtain.  Ve  believe  that  theory 
as  it  now  stands  is,  with  these  specific  restrictions,  applicable 
to  many  problems  of  more  than  passive  interest  to  DoD. 


4.  RADIATIVE  TRANSFER  IN  CONDUCTIVE  FIBER  CLOUDS 
4 .1  Basic  Theory. 


I 

The  purpose  of  the  radiative  transfer  portion  of  the  program 
is  two-fold.  First,  we  are  developing  a  simple,  efficient  computer 
program  to  handle  the  numerical  computation  of  radiative  transfer 
I  in  slabs  for  a  wide  range  of  s c a t t e r i ng / a b s o r b i ng  particle 

parameters,  and  for  quite  general  particles.  Second,  we  are 
integrating  the  RT  program  with  a  computer  program  that  treats 
the  scattering  and  absorption  specifically  by  lossy  conductive 
fibers,  in  order  to  determine  the  overall  electromagnetic  behavior 
of  such  configurations. 

The  basic  RT  theory  has  been  described  earlier.^®  Essen¬ 

tially,  we  are  using  the  van  de  Hulst  doubling  method  with  a  Taylor 
series  representation  for  the  starting  (thin-layer)  values  of 
the  transmission  and  reflection  matrices,  including  up  to  fourth- 
power  terms  in  the  optical  depth.  This  innovative  treatment 
'  of  starting  values  enables  us  to  avoid  the  cummulative  error 

build-up  incurred  when  one  begins  with  infinitesimal  starting 
layer  thicknesses. 

In  brief,  the  theory  runs  as  follows.  In  order  to  interpret 
the  transmission  and  reflection  matrices,  which  are  the  basic 
quantities  computed,  consider  Fig.  31.  The  layer  of  optical  depth 
z  contains  a  cloud  of  randomly  oriented  and  positioned  fibers  (or 
other  scattering  particles).  Radiat'ion  is  incident  from  outside 
the  layer,  and  we  characterize  the  right-going  and  left-going 
parts  at  any  point  in  the  layer  by  u  and  v,  respectively.  The 
spatial  direction  variable  has  been  discretized,  so  that  more 
specifically  u  is  an  N-component  column  vector  giving  the  angular 
distribution  of  radiation  into  the  right  hemisphere,  and  similarly 
for  v.  Now  u(0)  represents  radiation  incident  on  the  front  face 
of  the  layer,  giving  rise  to  transmitted  intensity 

u ( x )  =  T  u(0)  (43) 


and  reflected  intensity 


v  (  0 )  =  R  u  (  0 )  , 


(44) 


Figure  31.  Discretized  Right-Going  (u)  and 
Left-Going  (v)  Radiation. 

both  of  which  can  be  computed  once  the  transaiision  and  reflection 
Matrices  T  and  R  are  known. 


One  begins  by  evaluating  T  and  R  numerically  for  a  layer 
very  thin  compared  with  optical  depth,  i.e.,  x  <<  1,  using  the 

f  omul  a  s 

T  «  1  -  Qj  x  +  (1/2!)  (Qt2  +  Qz2)  x2 

-  (1/3!)  (Q^  +  Q2  Qj  Q2  +  2  Qj  Q22  +  2  Q22  Qt)  x3 

+  (1/4!)  {Q*  +  Q2  Qt2  Q2  +  3  Q22  +  3  Q22  Q  2 

+  3  Q1  °2  Q1  °2  +  3  °2  °1  °2  Q1  +  5  S'  +  5  Q,  S'  °i>  ^  (45) 

R  55  Q2  x  -  (1/2!)  (Qt  Qz  +  Qz  «  )  x2 

+  (1/3!)  (Qj2  Q2  +  Q2  Q 1 2  +  2  Q23  +  2  Ql  Q2  x3 

-  (1/4!)  (Qj3  Q2  +  Q2  Qa3  +  3  Qi  Q2  Ql2  +  3  Q  2  Q.,  Q 

+  3  Q  Q,  Q  2  +  3  Q  2  0,0^5  Q,  Q  3  +  5  Q  3  Q.)  x4  .  (46) 


These  equations  give  a  particularly  compact  method  for  obtaining 
T  and  R  to  high  accuracy)  the  scattering  behavior  of  the  fibers 
themselves  enters  in  the  determination  of  the  quantities  Q2  and 
Q2  •  which  are  defined  explicitly  in  terms  of  the  single  fiber 
scattering  pattern. 10 

Having  obtained  T  and  R  for  the  thin  layer,  the  corresponding 
values  for  layers  of  arbitrary  optical  depth  are  now  computed  by 
successive  applications  of  the  doubling  equations.  Vriting  T2  = 
T(2x),  R2  =  R(2x),  one  has 


T2  =  T  (1  -  R2)-l  T. 

R2  =  R  +  TR  (1  -  R2)-l  T  . 


(47) 

(48) 


There  are  thus  one  matrix  inversion  and  several  multiplications 
required  at  each  stage  of  doubling.  The  starting  layer  thickness 
should  be  chosen  judiciously  in  order  to  optimize  the  accuracy  of 
the  computation)  note  that  the  smaller  the  initial  choice  for  x, 
the  more  accurate  will  be  the  starting  values  for  T  and  R  from 
Eqs.  (45)  and  (46)  but.  on  the  other  hand,  the  more  repetitions 
of  the  doubling  equations  that  will  be  required,  with  attendant 
loss  of  precision,  in  order  to  attain  a  given  optical  depth. 
Numerical  trials  have  been  carried  out  which  ensure  that  we  work 
with  near-optimum  starting  depths. 

For  the  discretization  of  the  radiative  transfer  equation  we 
have  employed  the  extended  Simpson's  3/8  rule,  which  can  be 
derived  without  difficulty  from  the  ordinary  3/8  rule  available 
in  the  1 i t e r a tur e . 1 2  As  the  extended  rule  is  not  too  well  known, 
we  state  it  here: 


f  <izfU)  --  (3h/8)  [lU0)  +  3  fU^  +  3f(z2) 


+  2  +  3  )  i  •  •  •  + 1  (z.  _  )1  +  R , 


The  range  of  integration  ia  divided  into  3a  intervals  each  of 
width  h,  where  h  «  (a-b)/3s,  and  we  choose  a  odd,  so  that  3s  =  3, 
9,  15,  21,  ...  .  The  reaainder  tern  can  be  estimated  by 

R3s  =  -U4/80)  [{(Hi)  (b)  _  f(iii)  (t)j  (50) 

in  terns  of  the  third  derivatives  of  the  integrand  at  the  end 
points  . 

In  the  present  application,  a  is  identified  with  p  =  cos  G, 
where  G  is  the  angle  formed  by  the  radiation  flnz  with  the  forward 
direction  (z  azis  in  Fig.  31).  Thus,  T  and  R  will  determine 
radiation  patterns  over  the  discrete  set  of  directions  having 
uniform  increments  in  cos  G.  The  eztended  3/8  rale  enables  us  to 
avoid  the  equatorial  direction  G  =  n/2,  which  can  cause  numerical 
difficulties.  It  also  permits  us  to  readily  adjust  the  grid 
spacing  in  cos  G,  something  which  is  not  possible  with  the  more 
commonly  employed  Gaussian  quadrature,  because  of  the  complezity 
of  determining  weighting  coefficients  and  ordinates  in  the  latter. 


The  basic  single-particle  inputs  needed  to  the  RT  computer 

program  are  the  albedo 

scattering  cross-section 

eztinction  cross-section  *  '  ' 

and  the  phase  function  p0  (p,  p'),  obtained  by  simply  integrating 
the  scattered  power  p  (p,  0,  p',  f)  over  all  azimuths  0,  0*. 

Here  p  =  cos  0  and  0  are  the  angles  of  incidence,  and  p*  =  cos 
8',  0'  the  angles  of  observation. 

There  are  two  ways  of  proceeding,  depending  on  whether  or 
not  analytical  ezpressions  are  available  for  the  phase  function. 

The  two  simplest  analytical  cases  are!3 


po  ( P »  P  * )  =  1 


(Isotropic)  (S2) 


3  1 

Po<P*P'>  =  j  [1+P2p'2  +  J  (1-p2) (1-p2) (1-p'2) ] ,  (Rayleigh  Sphere) 


P0(p,p')  =  j  ll+p2p*2  +  j  ( 1-p2)  ( 1-p2  )  ( 1  —  |i  *2  )  ]  ,  (Rayleigh  Sphere) 


the  latter  being  appropriate  for  small  conducting  spheres  with 
unpolarized  incident  radiation.  These  cases  are  useful  for 
verifying  the  validity  and  accuracy  of  the  computer  program) 
numerical  results  are  readily  available  for  them  in  the  literature. 

In  the  absence  of  analytical  expressions  the  phase  function 
must  be  handled  in  tabular  form.  Following  Ishimaru,22  this 
is  done  as  follows.  For  inputs  we  employ  the  s  ingl  e- s  c  a  1 1  e  r  e  r 
intensities  S(pj)  evaluated  at  (3s  +  1)  equally-spaced  points 
in  p.  The  reflection  and  transmission  matrices  R  and  T  will 
each  have  dimensions  (3s  +  l)/2  x  (3s  +  l)/2.  Note  that  S(p) 
represents  the  scattered  intensity  averaged  over  all  scatterer 
orientations,  for  incidence  in  the  direction  p  =  pj  =  1  (the 
forward  direction). 

The  phase  function  can  now  be  determined  from 


P0(p,y ' ) 


2  Wm  Pm  (y)  P  ( y ’ ) 
m=0  m  m  m 


in  terms  of  the  Legendre  polynomials  PB.  The  latter  are  determined 
in  standard  fashion  using  a  recursion  sub-routine  based  on  the 


formulas 


p0  -  1 
Pi  -  H 


Pm+l(l,)  =  (stt)  [<2'n+1)lJ  V10-1"  pm-l<l,)] 


The  constants  Wm  are  given  by 


-  (-MC 


dy  P  (y)  S  (y)  , 


and  the  summation  of  Eq.  (54)  it  terminated  when  the  WB  become 
negligible.  The  integral  of  Eq.  (56)  is  approximated  using 
Simpson's  3/8  rule,  Eq.  (49),  with  grid  spacing  h  -  2/3s. 


The  doubling  method  can  now  be  applied  to  obtain  the  reflection 
and  transmission  matrices  for  layers  of  any  desired  optical 
length.  In  order  to  obtain  physical  results,  matrix  elements 
must  then  be  weighted  appropriately  to  allow  for  the  specific 
quadrature  rule  employed.  If  we  define  normalizing  constants 

c ( ti )  =  (an  Un)^,  n  =  1,2 . %  (3s  +  1)  ,  (57) 


where  the  a„  are  the  3/8-rule  weighting  coefficients,  then  the 
net  fluxes  reflected  in  the  back  hemisphere,  or  transmitted  in 
the  forward  hemisphere,  are  given  respectively  by 


%(3s+l) 

R  =  —prr-  y  R(n,  1)  c(n) 
c(1)n=l 

(5  b) 

%(3s+l) 

T  =  2  T(n ,  1 )  c(n)  . 


These  equations  are  for  radiation  incident  normally  on  the  slab, 
which  explains  why  only  the  first  columns  of  R  and  T  are  needed. 
For  isotropically  incident  radiation,  all  elements  of  R  and 
T  would  come  into  play.  For  the  angular  intensity  patterns,  on 
the  other  hand,  one  has 


R(n)  =  R(n , 1 ) / 2  c(l)c(n)  , 

T ( n )  =  T(n , 1 ) / 2  c(l)c(n)  ,  n  =  1,  2,  %  (3s+l) 


(59) 
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For  the  conductive  fibers,  the  first  step  in  the  analysis 
involves  orientstion  sversging  of  the  single-fiber  differential 
scattering  patterns.  If  p  (a,0)  is  a  unit  vector  along  the  fiber 
axis,  with  spherical  angles  a,0  relative  to  a  fixed  cartesian 
coordinate  system,  for  a  linearly  polarized  incident  wave  ejn 
along  the  positive  z  direction  we  have  o^iff  =  a^iff  (a, 0,9,0), 
and  the  or  ien ta t ion- ave r a ged  differential  scattering  pattern  is 
given  by 


’diff 


(2V  (  71 

(0)  =  (^)  J  daj  dBsinft  (ot,3,0,<f) 


Note  that  strictly  speaking  it  is  necessary  to  also  average  over 
the  polarization  angles  of  the  incident  and  scattered  waves. 
Because  we  are  only  interested  in  scalar  intensities,  however,  it 


incidence  (plane  formed  by  the  fiber  and  the  direction  of  inci¬ 
dence).  The  scattered  intensity  is  then  taken  without  regard  to 
its  polarization.  It  is  not  difficult  to  verify  that  this 
simplicatiou  does  not  affect  the  angular  dependence  of  the  result. 


The  scattering  patterns  are  actually  computed  in  fiber 
coordinates,  i.e.,  in  terms  of  the  angle  cos  6f  =  p*kOU{  formed 
by  the  fiber  and  the  observation  direction  kout  (9,0).  The  angle 
9f  can  be  expressed  in  terms  of  the  orientation  angles  by  noting 


that 


cos  9f  =  sin  9  sin  0  cos  (0-o)  +  cos  9  cos  0  . 


(61) 


We  see  that  the  0-dependence  drops  out  during  the  integration,  so 
that  the  resulting  pattern,  Eq.  (60),  is  rotationally  symmetric. 


The  short-fiber  limit  kh  <<  1  (h  =  half-length  of  fiber) 
provides  a  useful  check  on  the  computation.  In  this  limit  we  can 
assume  that  the  axial  current  induced  in  the  fiber  is  proportional 
to  ein*P  =  *in  0  and  that  the  corresponding  scattered  amplitude 
is  given  by  ein»p  sin  9f,  so  that 


°diff  =  °o  (ein*P  sin  9f)  *  (6^ 

where  a0  is  the  maximum  cross  section  obtained  at  broadside 
incidence  and  observation.  Pntting  this  expression  in  Eq.  (60). 
and  dividing  by  a  factor  of  4  to  correct  for  not  averaging  over 
incident  and  received  polar  ixa t ion,  gives  finally 

o,.„/oA  =  (1/30)  (3  +  cos20)  .  (63 
dill  o 

This  3  +  cos2  6  angular  dependence  constitutes  a  new  elemen¬ 
tary  scatterer  for  radiative  transfer,  which  might  be  called  the 
Rayleigh  fiber.  Note  that  it  falls  intermediate  relative  to 
the  isotropic  and  Rayleigh  sphere  (1  +  cos2  0)  cases. 

Numerical  computations  of  or ie n t a t i on- a v e r a g e d  scattering 
from  perfectly  conducting  fibers  are  shown  in  Fig.  32  for  several 
values  of  kh  (for  clarity  normalized  to  4  in  the  forward  direction 
P  =  cos  6  =  1).  The  uppermost  curve  is  the  Rayleigh  fiber 
pattern  of  Eq.  (63),  and  one  notes  that  the  numerical  results 
converge  smoothly  to  this  limit  for  kh  <<  1  (the  kh  »  0.1  curve, 
not  shown,  is  indistinguishable).  As  kh  becomes  greater  than 
unity  a  significant  peak  is  seen  to  develop  in  the  forward 
direction,  typical  behavior  for  scatterers  when  dimensions  are 
not  small  compared  to  wavelength. 

It  is  interesting  to  note  that  the  orientation-average  has 
been  computed  by  Borison,  but  for  half-wave  dipoles  ( kh  =  n/2).2* 
Upon  averaging  his  results  over  incident  and  received  polarizations 
one  obtains  precisely  Eq.  (63).  This  does  not  agree  with  the 
curve  of  Fig.  32»  we  believe  that  Borison' s  assumption  of  simple 
cosine  behavior  of  the  surface  current  for  all  fiber  orientations 
is  no  longer  adequate  at  kh  =  n/2. 


Now  consider  the  slab  geometry  for  a  cloud  of  randomly 
oriented  fibers,  as  shown  in  Fig.  33.  The  cloud  has  optical  depth 
x,  and  one  desires  to  compute  the  angular  intensities  transmitted 
and  reflected,  T(0)  and  R(9),  respectively,  as  shown  in  the 
figure. 

A  simple  check  is  available  on  numerical  computations  in  the 
thin-layer  limit  x  <<  1.  As  shown  in  Fig.  34,  the  transmitted 
intensity  per  steradian  in  the  direction  p  =  cos  0  must  originate 
in  the  differential  cone  of  fibers  indicated.  But,  in  the 
thin-layer  limit,  the  contribution  from  an  element  of  the  cone 
will  be  proportional  to  albedo  a,  the  phase  function  p(p),  and 
exponential  factors  accounting  for  attenuation  of  radiation  over 
the  travel  path  within  the  slab.  Thus,  with  x'  =  x-ps, 

fX/U  V-  e 

I(u)  =  Kup(u)  \  ds  e  e  ,  y  >  0  .  (64) 

Jo 

A  similar  equation  can  be  written  down  for  the  reflected  intensity, 
and  carrying  out  the  integrations  gives 

|,.-x  (1-e-<l-P)*/f)/(1.„)i  vt0 

I  (u)  =  >.'tp(u)  /  (65) 

(  (l-y)  ,  y  <  0 

These  last  formulas  are  given  by  van  de  Hulst*5  for  isotropic 

case  p(p)  =  1»  clearly,  from  our  derivation  they  must  continue  to 
hold  for  arbitrary  p(p). 

Intensity  patterns  have  been  computed  numerically  for  an 
optical  depth  x  =  0.1  and  the  three  elementary  scatterers, 

isotropic,  Rayleigh  sphere  (dipole),  and  Rayleigh  fiber.  Results 
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*re  shows  by  the  points  in  Fig.  35.  The  corresponding  snslyticsl 
approximations  of  Eq.  (65)  ere  given  by  the  solid  corves  in  the 
figure,  and  excellent  agreeaent  is  seen  throughout.  Incidentally, 
van  de  Hulst  gives  the  theoretical  value  <  =  1/4  for  the  propor¬ 
tionality  constant  of  Eqs.  (64)  and  (65).  Ve  find  that  auch  better 
results  are  obtained  using  a  somewhat  larger  value  obtained  by 
matching  Eq.  (65)  to  the  computed  point  for  noraal  reflection  for 
transmission  ( p  =  +  1).  For  exaaple,  for  isotropic  scatterers 
the  two  computations  then  agree  to  four  significant  figures, 
excepting  only  the  point  p  *  -.04762,  where  agreement  is  slightly 
poorer . 

4.5  Numerical  Results. 

A  plotter  routine  has  been  developed  for  the  reflected  and 
transmitted  intensity  patterns,  and  some  typical  results  will 
now  be  shown  for  normal  incidence  on  slabs  of  optical  thickness  x. 

In  Fig.  36  angular  patterns  are  shown  versus  observation 
direction  p  =  cos  6  for  lossless  isotropic  scatterers.  For 
this  case  s  =  9,  and  R  and  T  are  of  dimension  14  x  14.  Note 
that  our  quadrature  scheme  avoids  the  equatorial  plane  p  «  0 , 
where  analytical  difficulties  arise  (the  curves  are  easily  extra¬ 
polated  if  desired,  however).  In  the  transmitted  patterns  the 
coherent  field,  that  portion  of  the  originally  incident  beam 
which  has  survived,  is  represented  by  a  discrete  version  of 
the  delta  function  (sharply  rising  curves  near  p  =  1,  truncated 
at  unit  intensity  for  clarity).  As  was  observed  earlier**  these 
results,  including  the  discrete  delta  function  contribution, 
exactly  satisfy  energy  balance  constraints. 

The  same  computation  is  carried  out  in  Fig.  37,  but  for  much 
coarser  grid  spacing  (matrix  dimension  5  x  5).  Aside  from  the 
obvious  change  in  the  discrete  delta  function  contribution, 
these  curves  are  seen  to  be  effectively  identical  with  those 
of  Fig.  36.  This  demonstrates  the  stability  of  the  computer 
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program  with  respect  to  grid  spacing  at  least  for  this  simple 
case,  and  also  confirms  the  ease  with  which  grid  spacing  can 
he  adjusted  to  fit  the  problem  at  hand. 

In  Figs.  38  and  39  (albedo  »  0.8,  0.S  respectively)  losses 
are  introdnced,  with  corresponding  redaction  in  the  intensities. 
For  all  of  these  results,  checks  with  the  numerical  tables  of  van 
de  HulstlS  give  excellent  agreement. 

Figures  40,  41  and  42  show  the  intensity  distributions  obtained 
using  the  Rayleigh  phase  function  of  Eq.  (S3),  and  albedos  of 
1.0,  0.8,.  O.S.  Note  that  there  is  not  too  much  difference 
from  the  isotropic  case.  The  Rayleigh  case  differs  in  that 
single  psrticle  scattering  is  now  more  concentrated  in  the  forward 
and  back  directions.  In  terms  of  scattering  angle  0,  the  single¬ 
particle  intensity  is  given  by!3 

(1  +  cos^  0)  .  (66) 

the  effect  of  this  is  evident  when  one  compares  the  reflection 
curves,  which  have  become  relatively  somewhat  larger  in  the 
back  direction  (  p  =  -1),  although  not  much  change  has  occurred 
in  the  transmission  curves. 

In  Fig.  43  plots  of  the  transmitted  and  reflected  intensities 
are  given  for  an  optical  depth  of  unity,  where  the  single  scatterer 
is  1)  isotropic,  2)  Rayleigh  sphere,  3)  the  Rayleigh  fiber  of 
Eq.  (63),  and  4)  the  short  fiber  (parameters  as  noted  in  the 
figure).  Note  that,  just  as  was  true  for  the  phase  functions, 
the  radiative  transfer  intensities  for  the  Rayleigh  fiber  fall 
intermediate  between  the  isotropic  and  sphere  case.  It  is  also 
interesting  that  the  Rayleigh  fiber  and  short  dipole  (kh  =  0.5) 
results  are  now  indistinguishable,  even  though  their  phase 
functions  are  measurably  different  (see  Fig.  32). 
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Wavelength—  1.0E-06 
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kh-  5.0E-01 
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Dependence  of  the  ET  intensities  on  fiber  length  is  indicated 
in  Fig.  44,  for  fairly  highly  conducting  fibers.  From  the  figure 
one  can  see  that  both  the  transmitted  and  reflected  intensity, 
integrated  over  all  angles  as  in  Eq.(58),  show  a  mild  peak  at 
about  kh  =  3,  although  the  reason  for  this  is  not  clear. 

For  fixed  fiber  length  kh  =  3,  dependence  of  the  RT  inten¬ 
sities  on  optical  depth  x  is  shown  in  Fig.  45.  As  one  would 
expect,  the  reflected  patterns  increase  monotonically  with  x,  at 
the  same  time  becoming  more  nearly  isotropic.  The  transmitted 
patterns  show  the  same  behavior  in  the  diffusely  transmitted  flux 
up  to  about  x  =  1,  with  corresponding  reductions  in  the  coherent 
flux  (because  of  truncation  the  latter  effect  is  not  obvious  from 
the  figure).  With  further  increases  in  optical  depth  the  trans¬ 
mitted  flux  begins  to  decrease,  presumably  because  scattering  and 
absorption  mechanisms  now  dominate. 

Finally,  Fig.  46  illustrates  the  effect  of  varying  the  fiber 
conductivity,  maintaining  a  fixed  fiber  length  kh  *  1  and  optical 
depth  x  *  1 .  Both  reflected  and  transmitted  patterns  are  seen  to 
increase  monotonically  with  conductivity,  both  also  showing  signs 
of  saturation  for  the  uppermost  curves,  which  are  effectively 
approaching  the  perfectly  conducting  limit.  Again,  the  curves 
behave  qualitatively  as  would  be  expected,  in  view  of  the  fact 
that  the  albedo  is  increasing  toward  unity  with  increasing  conduc¬ 
tivity. 


The  above  examples  illustrate  the  usefulness  and  versatility 
of  the  computer  programs.  These  programs  are  now  sufficiently 
well  developed  to  begin  running  specific  problems  of  practical 
interest  to  CRDC  involving  both  single  fiber  scattering  and 
absorption,  and  the  radiative  transfer  properties  of  the  corre¬ 
sponding  fiber  clouds. 
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